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Abstract 

In this paper we present exact calculations of the partition function Z of the g-state Potts model and its 
generalization to real q, for arbitrary temperature on n- vertex strip graphs, of width L y = 2 and arbitrary 
length, of the triangular lattice with free, cyclic, and Mobius longitudinal boundary conditions. These parti- 
tion functions are equivalent to Tutte/Whitney polynomials for these graphs. The free energy is calculated 
exactly for the infinite-length limit of the graphs, and the thermodynamics is discussed. Considering the 
full generalization to arbitrary complex q and temperature, we determine the singular locus B in the corre- 
sponding C 2 space, arising as the accumulation set of partition function zeros as n — ► oo. In particular, we 
study the connection with the T = limit of the Potts antiferromagnet where B reduces to the accumulation 
set of chromatic zeros. Comparisons are made with our previous exact calculation of Potts model partition 
functions for the corresponding strips of the square lattice. Our present calculations yield, as special cases, 
several quantities of graph-theoretic interest. 
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1 Introduction 



The g-state Potts model has served as a valuable model for the study of phase transitions and critical 
phenomena fjl |j. On a lattice, or, more generally, on a (connected) graph G, at temperature T, this model 
is defined by the partition function 

W 

with the (zero-field) Hamiltonian 

n = -jJ2 6 ^ (i-2) 

(y) 

where <Ji = 1, ...,q are the spin variables on each vertex i G G; /3 = (fee? 1 ^ 1 ; and (jj) denotes pairs of 
adjacent vertices. The graph G = G(V, E) is defined by its vertex set V and its edge set E; we denote the 
number of vertices of G as n — n{G) = \V\ and the number of edges of G as e(G) = \E\. We use the notation 

K = /3J, a = vT 1 =e K , v = a-l (1.3) 

so that the physical ranges are (i) a > 1, i.e., u > corresponding to oo > T > for the Potts ferromagnet, 
and (ii) < a < 1, i.e., — 1 < v < 0, corresponding to < T < oo for the Potts antifcrromagnct. One defines 
the (reduced) free energy per site / = —j3F, where F is the actual free energy, via 

f({G},q,v) = lim ]n[Z(G,q t v)V n ] . (1.4) 

n — >oo 

where we use the symbol {G} to denote Hm^-nx, G for a given family of graphs. 

Let G 1 = (V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex set V and an edge 
set E' C E. Then Z(G, q, v) can be written as the sum 

Z(G, 9 ,«) = ^ g fc ( G ') w e(G') 
G'CG 

ti(G) e(G) 

= E £w« s (i.5) 

r=k(G) s=0 

where fc(G') denotes the number of connected components of G' and z rs > 0. Since we only consider 



connected graphs G, we have k{G) = 1. The formula (L5) enables one to generalize q from Z + to R+ 
(keeping u in its physical range). This generalization is sometimes denoted the random cluster model |^]; 
here we shall use the term "Potts model" to include both positive integral q as in the original formulation 



in eqs. (1.1) and (1.2), and the generalization to real (or complex) q, via eq. fll.q ). The formula (1.5) shows 
that Z(G, q, v) is a polynomial in q and v (equivalently, a) with maximum and minimum degrees indicated in 
cq. ( [l.5|) . The Potts model partition function on a graph G is essentially equivalent to the Tutte polynomial 
||-[pj and Whitney rank polynomial 0,0, (9|-@ for this graph, as discussed in the appendix. 



In this paper we shall present exact calculations of the Potts model partition function for strips of the 
triangular lattice of width L y = 2 vertices (or equivalently edges) and arbitrary length L x with various 
boundary conditions. This is a natural continuation of a previous study by one of us of the Potts model 
model on the analogous strips of the square lattice jl2[ ^3), and the reader is referred to that paper for 
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background and further references. We envision the strip of the triangular lattice as being formed by 
starting with a ladder graph, i.e. 2 x L x strip of the square lattice, and then adding edges joining, say, the 
lower left and upper right vertices of each square. The longitudinal (transverse) direction is taken as the 
horizontal, x (vertical, y) direction. We use free transverse boundary conditions and consider free, periodic 
(= cyclic), and Mobius longitudinal boundary conditions. These families of graphs are denoted, respectively, 
as S m (for open strip), L m (for ladder), and ML m (for Mobius ladder), where L x = m + 1 (edges) for S m 
(following our labelling convention in ]2(|) and L x — m for L m and ML m . One has n(S m ) = 2(m + 2) and 
n(L m ) = n{ML m ) = 2m. Each vertex on the cyclic strip L m has degree (coordination number) A = 4; this 
is also true of the interior vertices on the open strip S m . The Mobius strip can involve a seam, as discussed 
in Q. 

The motivations for these exact calculations of Potts model partition functions for strips of various 
lattices were discussed in Clearly, new exact calculations of Potts model partition functions are of 

value in their own right. In addition, it was shown Jl3| that these calculations can give insight into the 
complex-temperature phase diagram of the 2D Potts model on the given lattice. This is useful, since the 
2D Potts model has never been solved except in the q — 2 Ising case. Furthermore, with these exact results 
one can study both the T — critical point of the q-state Potts ferromagnet and, for certain q values (q = 2 
for the square strip; q — 2,3 for the present strip of the triangular lattice) the T = critical point of the 
Potts antiferromagnet. In addition, via the correspondence with the Tutte polynomial, our calculations yield 
several quantities of relevance to mathematical graph theory. 

Various special cases of the Potts model partition function are of interest. One special case is the zero- 
temperature limit of the Potts antiferromagnet (AF). For sufficiently large q, on a given lattice or graph 
G, this exhibits nonzero ground state entropy (without frustration). This is equivalent to a ground state 
degeneracy per site (vertex), W > 1, since So = InW. The T = (i.e., v = —1) partition function of the 
above-mentioned g-state Potts antiferromagnet (PAF) on G satisfies 

Z(G,q,-l) = P(G,q) (1.6) 

where P(G, q) is the chromatic polynomial (in q) expressing the number of ways of coloring the vertices of 
the graph G with q colors such that no two adjacent vertices have the same color || [l5l . The minimum 
number of colors necessary for this coloring is the chromatic number of G, denoted x(G). We have 

W({G},q)= lim PiG.q) 1 ^ (1.7) 



Using the formula (L5) for Z(G, q, v), one can generalize q from Z + not just to R+ but to C and a from 
its physical ferromagnetic and antiferromagnetic ranges 1 < a < oo and < a < 1 to a £ C. A subset of 
the zeros of Z in the two-complex dimensional space C 2 defined by the pair of variables (q, a) can form an 
accumulation set in the n — ► oo limit, denoted B, which is the continuous locus of points where the free 
energy is nonanalytic. This locus is determined as the solution to a certain {G}-dependent equation jl^, [l3|| . 
For a given value of a, one can consider this locus in the q plane, and we denote it as B q ({G}, a). In the 
special case a = (i.e., v = — 1) where the partition function is equal to the chromatic polynomial, the zeros 
in q are the chromatic zeros, and B q ({G},a — 0) is their continuous accumulation set in the n — * oo limit 
a series of papers we have given exact calculations of the chromatic polynomials and nonanalytic 
loci B q for various families of graphs (for further references on this a = special case, see |u|). With the 
exact Potts partition function for arbitrary temperature, one can study B q for a ^ and, for a given value 
of q, one can study the continuous accumulation set of the zeros of Z(G,q,v) in the a plane; this will be 
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denoted B a ({G} 7 q). It will often be convenient to consider the equivalent locus in the u — 1/a plane, namely 
B U ({G}, q). We shall sometimes write B q ({G}, a) simply as B q when {G} and a are clear from the context, 
and similarly with B a and B u . One gains a unified understanding of the separate loci B q ({G}) for fixed a 
and B a {{G}) for fixed q by relating these as different slices of the locus B in the C 2 space defined by (q, a) 
as we shall do here. 

Following the notation in JHJ and our other earlier works on B q ({G}) for a = 0, we denote the max- 
imal region in the complex q plane to which one can analytically continue the function W({G},q) from 
physical values where there is nonzero ground state entropy as R\ . The maximal value of q where B q 
intersects the (positive) real axis was labelled q c ({G}). Thus, region Ri includes the positive real axis for 
q > q c ({G}). Correspondingly, in our works on complex-temperature properties of spin models, we had 
labelled the complex-temperature extension (CTE) of the physical paramagnetic phase as (CTE)PM, which 
will simply be denoted PM here, the extension being understood, and similarly with ferromagnetic (FM) and 
antiferromagnetic (AFM); other complex-temperature phases, having no overlap with any physical phase, 
were denoted Oj (for "other"), with j indexing the particular phase |39 . Here we shall continue to use this 
notation for the respective slices of B in the q and a or u planes. 

We record some special values of Z(G, q, v) below, beginning with the q — special case 

Z(G,0,u) = (1.8) 

which implies that Z(G,q,v) has an overall factor of q. In general (and for all the graphs considered here), 
this is the only overall factor that it has. We also have 



G'CG 

For temperature T — oo, i.e., v = 0, 



Z(G, l,v)= ve[G ' ] = « e(G) • ( L9 ) 

G'CG 

Z(G,q,0)=q n ^ . (1.10) 



Further, 



Z(G,q,-l) = P(G,q) = 



X(G)-1 



IT 



U(G,q) (1.11) 



where U(G,q) is a polynomial in q of degree n(G) — x(G). 

We note some values of chromatic numbers for the strip graphs considered here: 

X(Srn) = 3 (1.12) 

lT \ f3 if m = mod 3 /-, 1Q s 

l. 4 it m = 1 or m = I mod 3 

and 

X (ML m ) = i. (1.14) 

Hence, for q — 1, 2, 3 we have 

Z(G,q,-l) = P(G,q) = for G = S m ,L m ,ML m and q = 1,2 (1.15) 

Z(G,3,-1) =P(G,3) = 3! for G = S m or G = L rn=0 rnod 3 (1.16) 
where G = L rn= o mo( i 3 means that this applies to L m for m = mod 3. For the graphs G = S m ,L m , 



and ML m , (i) the result fll.lq ) implies that for q = 2, Z(G,2,v) contains at least one power of the factor 
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(v + 1) = a; for q = 1, one already knows the form of Z(G, 1, f) from ( |l.9| ); (ii) the result ( 1.16| ) implies that 
for the cases where x(G) — 4, Z(G, 3, v) contains at least one power of (v + 1) as a factor. 



Another basic property, evident from eq. (1.5), is that (i) the zeros of Z(G, q, v) in q for real v and hence 
also the continuous accumulation set B q are invariant under the complex conjugation q — > q*; (ii) the zeros 
of Z(G, q, v) in v or equivalently a for real q and hence also the continuous accumulation set B a are invariant 
under the complex conjugation a — > a* . 

Just as the importance of noncommutative limits was shown in (eq. (1.9) of) |jl7f on chromatic polyno- 
mials, so also one encounters an analogous noncommutativity here for the general partition function (|l.5| ) of 
the Potts model for nonintegral q: at certain special points q s (typically q s — 0, l...,x(G)) one has 

lim lim Z(G, q, v) 1 ' 71 ^ lim lim Z(G, q, v) 1/n . (1.17) 

n — >oo q — >q s q — >q s n — >oo 



Because of this noncommutativity, the formal definition (1.4) is, in general, insufficient to define the free 
energy / at these special points q s ; it is necessary to specify the order of the limits that one uses in eq. 
(1.17). We denote the two definitions using different orders of limits as f qn and f nq : 



f nq {{G},q,v)= lim lim n' 1 In Z{G, q,v) (1.18) 

n — *oo q — >q 8 

f qn ({G},q,v) = lim lim rT 1 In Z(G, q, v) . (1.19) 

q — >q s n — ^oo 

In Ref. Jl7[ and our subsequent works on chromatic polynomials and the above-mentioned zero-temperature 
antiferromagnetic limit, it was convenient to use the ordering W({G},q s ) — \im cl ^ qs \im n ^ oc P(G,q) 1 ^ n 
since this avoids certain discontinuities in W that would be present with the opposite order of limits. In 
the present work on the full temperature-dependent Potts model partition function, we shall consider both 
orders of limits and comment on the differences where appropriate. Of course in discussions of the usual 



g-state Potts model (with positive integer q), one automatically uses the definition in eq. (1.1) with (|l.2| ) 
and no issue of orders of limits arises, as it does in the Potts model with real q. 



As a consequence of the noncommutativity (1.17), it follows that for the special set of points q = q s 



one must distinguish between (i) (B a ({G},q s )) nq , the continuous accumulation set of the zeros of Z(G,q,v) 
obtained by first setting q = q s and then taking n — > oo, and (ii) (B a ({G}, q s ))qn, the continuous accumulation 
set of the zeros of Z(G,q,v) obtained by first taking n — > oo, and then taking q — > q s . For these special 
points, 

(B a ({G}, q s )) nq + (B a ({G}, q s )) qn . (1.20) 



From eq. (1.8), it follows that for any G, 

exp(/„ g ) = for q = Q (1.21) 

and thus 

(B a ) nq = $ for q = 0. (1.22) 

However, for many families of graphs, including the circuit graph C n , and cyclic and Mobius strips of the 
square or triangular lattice, if we take n — oo first and then q — > 0, we find that (B u ) qn is nontrivial. 



Similarly, from (l.£) we have, for any G, 

(B a ) nq =% for q=l (1.23) 
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since all of the zeros of Z occur at the single discrete point a = (and in the case of a graph G with no edges, 
Z = 1 with no zeros). However, as the simple case of the circuit graph shows jL3|, (B u ) qn is, in general, 
nontrivial. 

As derived in |l3[ , a general form for the Potts model partition function for the strip graphs considered 
here, or more generally, for recursively defined families of graphs comprised of m repeated subunits (e.g. the 
columns of squares of height L y vertices that are repeated L x times to form an L x x L y strip of a regular 
lattice with some specified boundary conditions), is 

Z(G,g,«)=^c Gj (A Gj (< z , u )r (1.24) 
3=1 

where N\ depends on G. 

The Potts ferromagnet has a zero-temperature phase transition in the L x — > oo limit of the strip graphs 
considered here, and this has the consequence that for cyclic and Mobius boundary conditions, B passes 
through the T = point u = 0. It follows that B is noncompact in the a plane. Hence, it is usually more 
convenient to study the slice of B in the u = 1/a plane rather than the a plane. Since a — > oo as T — ► and 
Z diverges like a e ( G > in this limit, we shall use the reduced partition function Z r defined by 

Z r (G, q, v) = a- e{G) Z{G, q, v) = u< G *>Z(G, q, v) (1.25) 

which has the finite limit Z r — > 1 as T — > 0. For a general strip graph {G s ) m of type G s and length L x = m, 
we can write 

Z r ((G s ) m ,q,v) = a <: ^ , V; . ; iA f; A ( ; .,.„)'" (1.26) 

3=1 3=1 



with 



\ Gs , J ,u=u e « G ^/ m \ Gs , J . (1.27) 



For example, for the L y = 2 strips of the triangular lattice of interest here, with free transverse boundary 
conditions and either free or periodic longitudinal boundary conditions, we have e(S m ) = 4to + 5 and 
e(L m ) = 4m. (In the case of the cyclic strip with m = 2, L m degenerates in the sense that it contains 
multiple edges connecting some pairs of vertices, and similarly, for m = 1, L m contains both multiple edges 
and loops; the above formula applies to these cases if one takes care to count these multiple edges and loops, 
as is discussed further below.) 



2 Strip of Triangular Lattice with Free Longitudinal Boundary 
Conditions 

In this section we present the Potts model partition function Z(S m , q, v) for the L y = 2 strip of the triangular 
lattice S m with arbitrary length L x = m + 1 (i.e., containing m + 1 squares) and free transverse and 
longitudinal boundary conditions. One convenient way to express the results is in terms of a generating 
function: 

oo 

T(S,q,v,z)= Y / Z (Sm,q,v)z m ■ (2.1) 

rn—0 
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We have calculated this generating function using the deletion-contraction theorem for the corresponding 
Tutte polynomial T{S m , x, y) and then expressing the result in terms of the variables q and v. We find 

T ( S >9>"> Z )= V(S,q,v,z) (2 - 2) 

where 

M{S,q,v,z) =As,o + A s ,iz (2.3) 

with 

4s,0 = q(v 5 + 5w 4 + 8u 3 + 2v 3 q + I0v 2 q + 5vq 2 + q 3 ) (2.4) 

A S ,1 =~q(v + l) 2 (v + q) 3 v 2 (2.5) 

and 

V(S, q,v,z) = l- (v 4 + 4v 3 + 7v 2 + Aqv + q 2 )z + (v + l) 2 (v + qfv 2 z 2 . (2.6) 
(The generating function for the Tutte polynomial T(S m ,x,y) is given in the appendix.) Writing 

2 

V(S,q,V,z) = H(l-\sjz) (2.7) 

3=1 

we have 



A 



S,(l,2) 



Tsi2 ± (3w + v 2 + q)\/Rsi2 



(2.8) 



where 

T S12 = v 4 + Av 3 + 7v 2 + Aqv + q 2 (2.9) 

and 

Rsi2 = q 2 + 2qv - 2qv 2 + 5v 2 + 2v 3 + v 4 . (2.10) 

Ref. [p8f presented a formula to obtain the chromatic polynomial for a recursive family of graphs in the 
form of sums of powers of A^'s starting from the generating function, and the generalization of this to the 
full Potts model partition function was given in [^3|. Using this, we have 

7(Q „ „,\ _ ( A S,0>*S,1 + As,l) \m , (^S,oAg,2 + As,l) xm f)-\-\\ 

^\O m ,q,V) - — A SI H jt r r A S2 (^-J-J-J 

\*S,1 — *S,2) {AS,2 — AS,l) 

(which is symmetric under As,i <-» As, 2)- Although both the As,/s and the coefficient functions involve the 
square root \/R~si2 and are not polynomials in q and v, the theorem on symmetric functions of the roots 
of an algebraic equation pQ ] guarantees that Z(S rn , q, v) is a polynomial in q and v (as it must be by (1.5) 
since the coefficients of the powers of z in the equation ( |2.6| ) defining these Asj's are polynomials in these 
variables q and v. 

As will be shown below, for q > 2 (and for < q < 1, with the f qn and {B u ) qn definitions) the singular 
locus B u for this strip consists of arcs that do not separate the u plane into different regions, so that the 
PM phase and its complex-temperature extension occupy all of this plane, except for these arcs. For these 
ranges of q, the (reduced) free energy is given by 

/=^lnAg,i. (2.12) 



G 



This is analytic for all finite temperature, for both the ferromagnetic and antiferromagnetic sign of the spin- 
spin coupling J. The internal energy and specific heat can be calculated in a straightforward manner from 
the free energy (2.12); since the resultant expressions are somewhat cumbersome, we do not list them here. 



In contrast, for the range 1 < q < 2, the Potts antiferromagnet on the L x — > oo limit of this open strip does 
have a phase transition at finite temperature (see eq. ( 2.2.6| )), but this has unphysical features, including a 



specific heat and partition function Z that are negative for some range of temperature, and non-existence of 
a thermodynamic limit independent of boundary conditions. This will be discussed further below. 

We next consider the limits of zero temperature for the antiferromagnetic and ferromagnetic cases. Let 
us define 

*<*>=f^ =§<-<;>"-• <^> 

and P(Ck, q) is the chromatic polynomial for the circuit (cyclic) graph Ck with k vertices, 

P(C k ,q) = (q-l) k + (q-l)(-l) k (2.14) 
so that £>2 = lj Ds = q — 2, D4 = q 2 — 3q + 3, and so forth for other D^s. In the T = Potts antiferromagnet 



limit v = — 1, A51 = D 2 = (q — 2) 2 and As,a = 0, so that eq. (2.2) reduces to the generating function for 
the chromatic polynomial for this open strip 

Equivalently, the chromatic polynomial is 

P(S m ,q)=q(q-l)(q-2) 2 ^ . (2.16) 

For the ferromagnetic case with general q, in the low-temperature limit v — > 00, 

A s ,i = v 4 + 4v 3 + 0{v 2 ) , X s ,2 = v 2 + 2{q- l)v + 0(l) as v 00 (2.17) 

so that |As i| is never equal to IA52I m this limit, and hence B u does not pass through the origin of the u 
plane for the n — ► 00 limit of the open square strip: 

u = 0?B u ({S}). (2.18) 

In contrast, as will be shown below, B u does pass through u — for this strip with cyclic or Mobius boundary 
conditions. 

2.1 B q {{S}) for fixed a 

We start with the value a = corresponding to the Potts antiferromagnet at zero temperature. In this case, 
Z(S m , q,v = —1) = P(S m , (7), where this chromatic polynomial was given in eq. ( |2.16 ), and the continuous 



locus B q = 0. For a in the finite-temperature antiferromagnetic range < a < 1, B q consists of a single 
self-conjugate arc crossing the positive real q axis at 

q c ({S}) = -v(3 + v) = (1 - o)(2 + a) (2.1.1) 



where the factor (3v+v +q) multiplying the square root in eq. (2.8) vanishes, so that the two roots coincide 



Parenthetically, we observe that this is same as the value of q c for the Potts model on the n — > 00 limit of 
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the circuit graph, C n [£4 fL3|. The arc endpoints occur at the branch points where the square root is zero, 



qs,end P t. = (a - l)[o - 2 ± 2ija } . (2.1.2) 

As a increases from to 1, the crossing point (2.1.1) decreases monotonically from 2 to 0, and as a reaches 
the infinite-temperature value 1, B q shrinks to a point at the origin. In the ferromagnetic range a > 1 
the self-conjugate arc crosses the negative real axis, since q c ({S}) < 0, and has support in the Re(q) < 
half-plane. In Figs. |l| and ^ we show B q and associated zeros of Z in the q plane for the antiferromagnetic 
value a = 0.1 and the ferromagnetic value a = 2. 



Im(q) 




-1 1 ' ' ' ' ' ' ' 1 

0.5 1 1.5 2 

Re(q) 

Figure 1: Locus Bq for the n — > 00 limit of the L y = 2 triangular strip {S} with free longitudinal boundary conditions, for 
a = 0.1. Zeros of Z(S m ,q, v = —0.9) for m = 20 (i.e., n = 44 vertices, so that Z is a polynomial of degree 44 in q) are shown 
for comparison. 



2.2 B u for fixed q 

We show several plots of the locus B u for various values of q in Figs. |^ - ||. For values of q where non- 
commutativity occurs, we consider B qn . Given the algebraic structure of A<jj, j — 1,2, the degeneracy 
of magnitudes |As,i| = \Xs,2\ and hence the locus B u occurs where (i) Tsi2 = 0, (ii) the prefactor of the 
square root vanishes: 3v + v 2 + q = 0, (iii) Rsu = 0, (iv) if q is real, where Rsi2 < so that the square 
root is pure imaginary, and (v) elsewhere for complex u, where the degeneracy condition is satisfied. All 
of these five possibilities are realized here, although some, such as (i) can yield solutions already subsumed 
by other conditions; for example, for the case q = 1.8, shown in Fig. ^, T$i2 vanishes at two real points, 
u ~ — 29.1, —1.35, both of which are subsumed within the solution of condition (iv), which is a line segment 
—69.486 < u < —1.291, and so forth for various other solutions. In cases where B M does not enclose regions, 
Ag.i is dominant everywhere in the u plane, and degenerate in magnitude with As ; 2 on B u ; the cases where 
B u does enclose regions will be discussed individually. 
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4 



lm(q) 



-3 -2 

Re(q) 



Figure 2: Locus B q : same as Fig. hi for a = 2 (i.e., v = 1). 



For large values of q, we find that S u consists of the union of (i) a circular arc centered at u ~ — 1/g that 
crosses the negative real u axis, and (ii) a line segment on the negative u axis. For example, for q = 10, the 
arc has its endpoints at two of the branch point zeros of y/Rsi2, at u = (7± >/l5i)/32 ~ 0.21875 ± 0.12103i, 
and the line segment, which is the solution to the condition (iv) above, that Rsi2 < for real q, has left and 
right endpoints at iti = — 1 and u% — —1/4. As q decreases, u\ and ui move toward more negative values, 
and as q decreases to 2, u\ — ► — oo, i.e. a — 0, while u-i — ► — 1. For q = 2, S„ is the union of the above line 
segment — oo < u < — 1 and the circular arc 



= -1 + V2e ±z 



< 9 <TT 



where 



arctan 



(^) - 13.5° . 



V 11 J 



That is, the endpoints of the circular arc occur at 



(2.2.1) 
(2.2.2) 

(2.2.3) 



Proceeding, we observe that for 1 < q < 2, the locus B u separates the u plane into different regions, while 
for < q < 1, as was true for the range q > 2, B u does not separate the u plane into different regions. In 
the range q m < q < 2, where 

'5\ 2 



Q = 1-n 



(?)■ 



1.5625 



(2.2.4) 



B u crosses the real u axis vertically at three different points and separates the q plane into three regions: 

(i) the paramagnetic phase, which includes the infinite-temperature point u = 1, where Xs.i is dominant; 

(ii) the surrounding O phase that extends to the circle at infinity in the u plane, i.e. to the point a = 0, 
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in which phase As,2 is dominant, and (iii) a third phase, of O type, located slightly to the left of u = 0, in 
which As. 2 is dominant. It is interesting that although As,2 is dominant in both of the O phases, there is 
still a boundary that separates them completely; this is a result of the fact that on this boundary, the other 
A, namely A5.1, becomes degenerate in magnitude with As. 2- Two of the points where B u crosses the real 
axis occur where the condition (ii) holds, i.e., where the prefactor 3v + v 2 + q multiplying the square root in 



(2.2.5) 



eq. (^81), vanishes, at v = (l/2)[-3± y/9 - Aq ], i.e., 

1 ± ^9 - 4 g 



2(2-g) 



For q = 2 — e with < e << 1, these crossings are at u = — 1 + O(e) and u ~ 1/e. As q decreases, the 
smaller (larger) solution moves to the right (left). Specifically, as q decreases from 2 to 0, the larger solution 
decreases monotonically from infinity to 1 while the smaller one increases monotonically from —1 to —1/2. 



For example, for q = 1.8, the crossings in eq. (|2.2.5[) are at u — (l/2)(5 ± 3a/5 ) ~ -0.854, 5.854. The 



existence of the right crossing on the positive real u axis for u > 1 means that the free energy of the Potts 
antiferromagnet is nonanalytic at the temperature 

T p .s = f -, for < q < 2 (2.2.6) 



kn hi 



(l/2)(-l + v /9^) 



(where both J and the log are negative). However, just as found in p3[ , this nonanalyticity has associated 
unphysical features, including negative Z , negative specific heat in the low-temperature phase, and non- 
existence of a thermodynamic limit independent of boundary conditions in the low-temperature phase. As 
q decreases from 2 to the value q rn given above, the left and right endpoints of the line segment merge, 



and it shrinks to a point at u\ — U2 — —4. For this value q — q m , the crossings in eq. ( 2.2.5 ) are 
u = (4/7)(2 ± \/lT ) — 3.038, —0.7524. In Fig. || we show B u for q = q rn . The arc endpoints occur at 
u = (20 ± 8V6* )/49 ~ 0.4082 ± 0.3999i. For q m < q < 2, there is a multiple point on the negative real axis 
where a branch of B u crosses this axis vertically and intersects with the line segment. There are also multiple 
points at complex-conjugate values of u where the circular arc intersects the closed curve. As q decreases 
from 2 to q mi the circular arc in the left-hand complex plane enlarges while the closed curve extending into 
the right-hand plane shrinks and becomes convex. As q decreases below q m , the closed curve evident to 
the left in Fig. ^ breaks apart, forming two complex-conjugate arcs, and in the interval 1 < q < q m , B u 
no longer contains any line segment but instead consists of these complex-conjugate arcs and the closed 
curve to the right; an illustrative example is given in Fig. [?], where the arc endpoints occur at u — ±2i and 
u = (l/9)(4± 2y/5i ) ~ 0.444 ± 0.497i. As q decreases toward 1, the arcs shrink to points at u — e ±47T / 3 . 
At q = 1, {B u ) qn is an oval (the solution to the degeneracy equation |u(l — u)\ = 1) that (i) crosses the real 
axis at u = (1/2)(1 ± \/5 ), i.e. at approximately, 1.618 and —0.618, and (ii) crosses the imaginary axis at 
u = ±[(l/2)(— 1 + \/Z )] 1 ^ 2 i — ±0.7862i. For < q < 1, B consists only of two disjoint arcs, as illustrated 
for the value q = 1/2 in Fig. ||. For q = 0, the locus (B u ) qn is the circular arc 

u=\e ie , ^<e<^ (2.2.7) 
2 2 ~ ~ 2 y J 

that crosses the real axis at u = — 1/2 and has endpoints at u = ±i/2. This is qualitatively similar to (B u ) qn 
for the L y = 2 open strip of the square lattice at q = 0, which was another circular arc crossing the negative 
real axis at u = —1/3 with endpoints at u = (— 1 ± 2\/2i)/9 H . 
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Note that at q = 0, 1, one encounters the noncommutativity ( 1.17 ); if one sets q to either of these values 
first and then takes n — > oo, the resultant B u = 0. 

Our findings for this L y = 2 open strip of the triangular lattice may be contrasted with those for the 
L y = 2 open strip of the square lattice in |H| . In the latter case, B u consisted of arcs that never enclosed 
any regions for real q > 0. 



0.4 




lm(u) o 




-0.4 1 1 1 1 1 1 1 ■ 1 1 1 1 1 1 1 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 

Re(u) 



Figure 3: Locus B u for the n — > oo limit of the L y = 2 triangular strip, with free longitudinal boundary conditions, {S} with 
q = 10. Zeros of Z(S m , q = 10, v) in u for m = 20 are shown for comparison. 
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Im(u) 



-6 -4 -2 2 4 

Re(u) 



Figure 5: Locus B u - same as in Fig. U for q = 1. 



lm(u) o 



V 



-4-3-2-10 1 2 3 4 

Re(u) 

Fi gure 6; Locus 3 U '- same as in Fig. ^] for q = (5/4)^. 
2 



lm(u) o - 




Re(u) 



Figure 7: Locus £>„: same as in Fig. H for q = 1.25. 
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-1 -0.5 0.5 1 1.5 2 

Re(u) 

Figure 8: Locus B u - same as in Fig. ^ for q = 0.5. 

3 Cyclic and Mobius Strips of the Triangular Lattice 

3.1 Results for Z 

By either using an iterative application of the deletion-contraction theorem for Tutte polynomials and con- 
verting the result to Z , or by using a transfer matrix method (in which one starts with a q 2 x q 2 transfer 
matrix and generalizes to arbitrary q), one can calculate the partition function for the cyclic and Mobius 
ladder graphs of arbitrary length, Z(G, q,v), G = L m , ML m . We have used both methods as checks on the 



calculation. Our results have the general form ( 1.24 ) with N\ = 6: 



Z(G m ,q, v) = J2 c GA x G,j(q,v)) m ,G m = L m , ML m (3.1.1) 

where 

Alj = A M z,,j ,j = 1,...,6 (3.1.2) 

We find 

Xl,i=v 2 (3.1.3) 

and 

Al,5 = As,i , Al : 6 = A5,2 . (3.1.4) 
The Xj for j = 2, 3, 4, are the solutions of the cubic equation 

r] 3 - v{v 3 + 8v + Av 2 + 2q)r] 2 + v 2 (2v 3 q + 6v 3 + 8qv + q 2 + 8v 2 + 2v 4 + 6qv 2 )r] 

- v 4 {v + l) 2 {v + q) 2 = . (3.1.5) 
For the cyclic strip the coefficients in eq. ( ft.l.lj ) are 

c L ,i = q 2 - 3q + 1 (3.1.6) 

c LJ = q-l for j = 2,3,4 (3.1.7) 

c LJ = l for 3=5,6 (3.1.8) 
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For the Mobius strip, we have 

cmlj = -1 (3.1.9) 

and 

c M lj = 1 for j = 5,6 . (3.1.10) 

The coefficients cmlj, j — 2,3,4 are more complicated and are given by the generating function (see the 
appendix) by the formula given in pa] . 



Our exact calculations (cf. eq. ( |3 . 1 . 2| ) ) yield the following general result 

B({L})=B({ML}) . (3.1.11) 

This is the same result that one of us found for the analogous strip of the square lattice [|l2| [ll| and is in 
accord with the conclusion that the singular locus is the same for an infinite-length finite-width strip graph 
for given transverse boundary conditions, independent of the longitudinal boundary condition. Owing to the 
equality ( fo.l.ll ), we shall henceforth, for brevity of notation, refer to both B({L}) and B({ML}) as B({L}) 



and similarly for specific points on B, such as q c ({L}) = q c ({M L}) , etc. 

Our main interest here is in large m and the m — > oo limit. However, for completeness, we make the 
following remark. If m > 3, then L m is a (proper) graph, but the m = 1 and m = 2 cases requires special 
consideration; in these cases, L m degenerates and is not a proper graph[]. L2 is the multigraph obtained 
from the complete graph K4 ^ by doubling two non-adjacent edges (i.e., edges that do not connect to any 
common vertex). L\ is the pseudograph obtained by connecting two vertices with a double edge and adding 
a loop to each vertex. Our calculation of Z(L m ,q,v) and the corresponding Tutte polynomial T(L m ,x,y) 
apply not just for the uniform cases m > 3 but also for the special cases m = 1, 2 if for m — 2 one includes 



the multiple edges and for m = 1 the multiple edges and loops in the evaluation of (1.1), (1.2), and (|1.5|). 
Note that in the T = case for the antiferromagnet, the resulting partition function, or equivalently, the 
chromatic polynomial, is not sensitive to multiple edges, i.e. is the same for a graph in which two vertices 
are connected by one edge or multiple edges; however, the general partition function (Tutte polynomial) is 
sensitive to multiple edges. The chromatic polynomial is sensitive to loops and vanishes identically when a 
pseudograph has any loops. 

3.2 Special values and expansions of A's 

We discuss some special cases. First, for the zero-temperature Potts antiferromagnet, i.e. the case a = 
(v = —1), the partition functions Z(L m , q, v) and Z(ML m , q, v) reduce, in accordance with the general result 



(1.6), to the respective chromatic polynomials P(L m , q) and P{ML m , q) |32J with Aj's comprised of the four 
terms 1, (q — 2) 2 , and (l/2)[5 — 2q ± y/9 — 4q\. The two remaining Al^'s vanish in this limit. (Since only 
four Al.j's occur in the chromatic polynomial, a different numbering convention was used in p2| than here, 
where, in general, six occur.) For the infinite-temperature value a = 1, we have Xlj — for j = 1, 2, 3, 4, 6, 
while Al_5 = q 2 , so that Z(G,q,a = 1) = q 2m = q n for G — L m , ML m , in accord with the general result 



(1.10). For the real interval q > 3 and the region R\ to which one can analytically continue from this interval 



(see Fig. || below), W = q — 2. Hence, W — 1 at q — 3. A technical remark is the following: one can, and 



it is convenient to, take the n — > 00 limit with m = mod 3, so that, by eqs. (1.13) and (fLTq), P = 3!, so 



1 A proper graph has no multiple edges or loops, where a loop is an edge that connects a vertex to itself. A multigraph may 
contain multiple edges, but no loops, while a pseudograph may contain both multiple edges and loops . 
2 The complete graph K p is the graph with p vertices each of which is adjacent to all of the other vertices. 
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that the limit for W exists at a — as well as at a ^ 0. In contrast, if one took n — ► oo using all positive 
integer values of m, then, strictly speaking, the limit for W in eq. (1.7) would not exist for a — 0, since 
P would have the nonconvergent values 6, 0, 0, 6, , .. for m = 3k, 3k + 1, 3fc + 2, 3k + 3, .. For the Mobius 
longitudinal boundary condition, no such convenient choice is possible, since x = 4 for all m; here there is 
a noncommutativity: if one starts with q slightly larger than 3, takes n — > oo to calculate W, and then lets 
q\3, one gets W(q — 3) = 1, but if one sets q — 3 first and then lets n — ► oo, one gets W(q = 3) = 0. 
At q = (with appropriate choices of branch cuts) we find that 



A 



LA 



A 



L,3 



Al : 2 — — — 



L,5 



v 2 + Av + 7 + (v + 3)\Jv 2 + 2v + 5 



and 



v 2 + Av + 7 - (v + 3)vV + 2v + 5 



Since there are dominant terms that are degenerate, namely A^ 2 = Aj^g, it follows that 

g = is on 6 g ({i}) Vs. 



(3.2.1) 
(3.2.2) 

(3.2.3) 
(3.2.4) 



This was also true of the circuit graph and cyclic and Mobius square strips with L y = 2 for which the 
general Potts model partition function (Tutte polynomial) was calculated in |l3|]. For q — 0, the coefficients 

1 for j = 2, 3, 4 so that the equal terms cancel each other pairwise, yielding 

The noncommutativity ( 1.17 ) occurs here: 



Cj = 1 for j = 1, 5, 6 and Cj 

Z{L m , q = 0, v) = 0, in accordance with the general result (1 
exp(/ ng ) = 0, while |exp(/ ?n )| = \\l,5\ 1/2 - 

At q = 1,2 we again encounter noncommutativity in the calculation of the free energy. For q = 1, 
f nq — 2 lna = 2K, while f qn depends on which phase one is in for a given value of a. For < a < a c (q = 1), 
where, from eq. J3.3.4Q , a c (q = 1) = (l/2)(-l + VE) ~ 0.6180, /,„ = (l/2)lnA L , c , where A L , C is the 
cube root that is dominant in this phase (corresponding to (1/2) [5 — 2q + -y/9 — 4gJ for a — > 0), while for 
a > a c (q = 1), f qn = (l/2)lnA L ,5- 

Similarly, for q = 2, again with an appropriate choice of branch cuts, 



A 



£,(2,4) 



v(v + 1) 



Al.i = A L: 3 = v 2 

! + 3v + 4 ± + l)(t; 2 + 5v + 8) 



1/2 



and 



A 



£,(5,6) 



(w + l)(v + 2) 



(u 2 +w + 2)± (w + l)(u + 2)(v 2 -u + 2) 



1/2 



(3.2.5) 
(3.2.6) 

(3.2.7) 



For this value, q = 2, the coefficients are c^^ = —1, while clj = 1, 2 < j < 6; hence, the (Al,i)™ and 
(A_L,3) m terms cancel each other and make no contribution to Z, which reduces to 



Z(L m ,q=2,v)= Yl (^j)' 

j=2,4,5,6 



(3.2.8) 



Hence also, f qn ^ f nq at g = 2. 

We observe that the Aj 's have a more symmetric form for q = 2 when expressed in terms of the variable 

u: 

\lxu = \l^u = u 2 {1-u) 2 (3.2.9) 
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A 



L,(2,A),u 



and 



A 



L,(5,6),« 



(1-U) 



(1 + u) 



1-u- 



2u 2 ± 



2u 2 ± 



(1 -u)(l + 3it + 4u 2 ) 



(l + u)(l-3u + 4u 2 ) 



1/2 



1/2 



(3.2.10) 



(3.2.11) 



One sees that each member of the pair A^ 



(2,4),« 



is equal to the respective member of the pair X L 



(5,6), 



with 



the replacement u — > — u. It follows that |Al,2,«| = IAl^iJ and |Al,4,«| = | Ai, 3 e,i*| f° r u pure imaginary. 

In general, for q G Z + , the partition function Z(L m ,q,v) for the cyclic width L y = 2 strip of the 
triangular lattice is identical to the partition function for the ID g-state Potts model with nearest-neighbor 
and next-nearest-neighbor spin-spin couplings that are equal in magnitude. This equality can be seen easily 
by redrawing the strip of the triangular lattice as a line with additional couplings between next-nearest- 
neighbor vertices on this line. In Ref. |4^] Tsai and one of us calculated Z for the latter model (with, in 
general, unequal nearest and next-nearest-neighbor spin-spin couplings). Hence, in particular, the q — 2 and 
q = 3 of eq. ( 3.1.1 ) coincide with the results in (section IX of) |f42f . 



In order to study the zero-temperature critical point in the ferromagnetic case and also the properties 
of the complex-temperature phase diagram, we calculate the Agj>'s corresponding to the Xaj's, using eq. 
(1.27). In the vicinity of the point u = we have 



A 



L.l.u 



= u 2 (l-uf 



and the Taylor series expansions 



Xl,2,u - 1 - 2u 3 + 2(q - l)u 4 + 0(u 5 ) 
Xl,3,u = u 2 + 0{u 3 ) 
Xl,4,u = u 2 + 0{u 3 ) 

Al, 5 ,„ = l + 2(q- l)u 3 l + u + 0{u 2 ) 
Al,6,« = u 2 + 2(q - 2)u 3 + 0(u 4 ) . 



(3.2.12) 



(3.2.13) 

(3.2.14) 
(3.2.15) 
(3.2.16) 
(3.2.17) 



Hence, at u — 0, \l,2,u an d Al,5,u are dominant and |Az,,2,u| = |Al,5,u|, so that the point u = is on B u for 



any q ^ 0, 1, where the noncommutativity (1.17) occurs. To determine the angles at which the branches of 
B u cross each other at u = 0, we write u in polar coordinates as u — re lS , expand the degeneracy equation 
|Al,2.m| = |Al.5,m|, for small r, and obtain qr 3 cos(3#) = 0, which implies that (for g ^ 0, 1) in the limit as 
r = \u\ 0, 

a _ (2j + IK 



6 



J = 0,1,..., 5 



(3.2.18) 



or equivalently, 9 — ±7r/6, ±7r/2, and = ±57r/6. Hence there are six curves forming three branches of B u 
intersecting at u — and successive branches cross at an angle of ir/3. The point u = is thus a multiple 
point on the algebraic curve B u , in the technical terminology of algebraic geometry (i.e., a point where 
several branches of an algebraic curve cross [0). In the vicinity of the origin, u = 0, these branches define 
six complex-temperature phases: the paramagnetic (PM) phase for — 7r/6 < 9 < 7r/6, together with the 
phases Oj for 1 < j < 5, with Oj occupying the sector (2j — 1)tv/6 < 8 < (2j + l)n/6. Note that O3 =03, 
O4 = O2, and O5 = 0;j\ For the case of interest here, namely, q > 0, \l,5,u is dominant in the PM phase 
and in the O2 and 0\ phases, while Al,2,« is dominant in the 0\, O^, and O3 phases. 
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One of our interesting findings is that for q = 2 and for q = 3 the Potts antiferromagnet on the infinite- 
length, width L y = 2 strip of the triangular lattice has a zero-temperature critical point. (As must be the 
case for this to be physical, this is independent of the longitudinal boundary conditions.) In the q = 2 Ising 
case, this involves frustration, and the 2D Ising antiferromagnet on the triangular lattice also has a T = 
critical point Q]. In contrast, for the q = 3 case, the T = critical point does not involve any frustration, 
and the q value at which this occurs for our L y = 2 strip is different than the value, q = 4, where the Potts 
antiferromagnet has a T = critical point on the full 2D triangular lattice. In order to study the T = 
critical point for the L y = 2 strip for these two values of q, it is useful to calculate expansions of the Xj's; only 
Xl,5 and Al,2 are necessary for physical thermodynamic properties, while the full set of 3 — lj 2, ■■, 6 is 
necessary for the study of the singular locus B. 



For q = 2, besides the exact expressions Xl,i = Ax 3 = (1 — a) 2 from eq. (3.2.5), we have the expansions 



A 



£,(2,4) 



-(1 



± iVa -a + -a 2 + 0(a 3 ) 



and 



a 2 

Al,(5,6) = a + y (1 + a 2 ) ± x/a 



0(a 3 



(3.2.19) 



(3.2.20) 



Note that (i) these are not Taylor series expansions in a, but rather in the variable y/a, and (ii) each member 
of the pair (2,4) is equal to the respective member of the pair A^ ( 5 6 ) with the replacement a — ► —a. 
As shown above, for f nq and B nq , where one sets q = 2 first and then takes n — ► 00, Alj, j = 1,3, make 
no contribution, and S ng is determined by the degeneracy in magnitude of Xlj, j = 2,4,5,6. From the 



expansions ( |3 .2.1 9| ) and (3.2.2C), it follows that in the neighborhood of the point a = 0, (B a ) nq is determined 



by the equation |Al,2| = \Xl,s\ and is a vertical line segment along the imaginary a axis. This will be 
discussed further below (see Fig. pHI) . 

For q = 3, besides the exact expression Al.i = (1 — a) 2 , we calculate the expansions 



Xl,(2,3) — e 



2e ±"/3 a + 

(where the upper (lower) sign applies for j = 2 (j = 3)); 



-1± 



2iV3 



)a 2 + 0(a 3 ) 



and 



A i>4 = 4a 2 - 4a 3 + 0(a 4 ) 
Al,5 = 1 + 6a - 3a 2 + (9(a 3 



Al, 6 = 4a 2 - 28a 3 + 0(a 4 



There are thus four dominant Alj's at a = 0, viz., those with j = 1, 2, 3, 5, which satisfy |Al ; 
a = pe 1 ^ and expanding the dominant Alj's in the neighborhood of a = 0, we obtain 

|Al,i| 2 = 1 -4p cos (j) + 0(p 2 ) 

\Xl,(2. 3 )\ 2 = 1 - 4pcos(0 T tt/3) + 0(p 2 ) 

and 

|A il5 | 2 = 1 + I2pcos0 + O(p 2 ) . 



(3.2.21) 



(3.2.22) 
(3.2.23) 

(3.2.24) 
1. Writing 



(3.2.25) 
(3.2.26) 

(3.2.27) 
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Equating dominant |Alj|'s, it follows that in the neighborhood of a — 0, B a consists of four curves, forming 
complex-conjugate pairs, passing through a = at the angles 



cf> = ±arctan 

and 



(^p) ~ 76.10° . (3.2.28) 



In Fig. H below we show B a for this case. 



= ±y = ±150° (3.2.29) 



3.3 B q ({L}) for fixed a 

3.3.1 Antiferromagnetic case < a < 1 

We begin with the case a = 0, i.e. the T = limit of the Potts antiferromagnet. The locus B q for this case 
was stated in [31 and given as Fig. 3 in j3^[. This locus B q , shown as Fig. [| here, separates the complex 
q plane into three different regions: (i) including the segments q > 3 and q < on the real axis, where 
W = (q — 2); (ii) R 2 , including the real interval 2 < q < 3, in which |W| = 1; and (hi) i?3, including the real 
segment < q < 2, in which \W\ — [(l/2)(5 — 2q + y/9 — 4q )] 1 / 2 . At q — q c there are actually four terms 
that are degenerate in magnitude, 

Al,i = A L , 5 = |A L , 2 | = |A Ll3 | = 1 at q = 3 for a = (3.3.1) 

corresponding to the property that this is a multiple point on B q (in the terminology of algebraic geometry) 
where four curves intersect. In contrast, the other two points at which B q crosses the real axis, at q = and 
q = 2, involve only the minimal number (two) of degenerate magnitudes and are hence not multiple points. 
Evidently, q c ({L}) = 3 for this a = case. 

We show calculations of B q ({L}) in Figs. [l(], [ll], and |l2| for finite temperature. As a increases from 0, 
rather than having a fourfold degeneracy of A's at q = 3, as in the a — case, eq. (3.3.1), one only has a 
twofold degeneracy, 

|Ai, 5 | = [Ai,i| (3.3.2) 



This occurs at the point 



, r in ,r fl-a)(3 + 2o) v(2v + 5) .„ „ „, 

q c ({L}) - q c ({ML}) = ± >- = - V ■ J . (3.3.3) 



1 + a 

Corresponding to this value of q is the pair 

1 



a c ({L})± = a c ({ML}) ± = - 



(q+l)± V? 2 -6q + 25 



(3.3.4) 



As a increases from to 1, q c ({L}) decreases monotonically from 3 to 0. Over the interval < a < ad, where 

a d = i(-3 + VT7 ) = 0.280776.. (3.3.5) 

the innermost region, i?2 continues to exist, but shrinks in size. Its left-hand boundary on the real q axis, 
where it is contiguous with region i?3, continues to lie at 

QR2-R3 = 2 f° r < a < a d ■ (3.3.6) 
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This follows from the fact that this point is the real solution to the degeneracy equation of leading terms 
|-Vl,i| = |<VmIj and these are both equal (to v 2 ) at q — 2. As a increases through the value a — ad, 
q c {{L}) a =a d — 2, so this innermost region R2 disappears; this can be seen from the fact that its right-hand 
boundary point q c becomes equal to its left-hand boundary point, given by qn 2 -R a . Thus, for a<j < a < 1, 
the locus B q separates the q plane into only two regions, R\ and R3. 

The Potts antiferromagnet has a phase transition at the temperature given by the upper choice of the 



sign in eq. ( 3.3.4 ), i.e., 

T P- L = I 1 ( J (fTX\ S ' ° < q < 3 ( 3 - 37 ) 

k B ln(a c ({L})+) 

(where both J and the log are negative). However, just as was found for the analogous phase transition 
of the Potts model on the infinite-length L y = 2 strip of the square lattice (which occurred for the range 
< q < 2 except for f nq for the integral value q = 1) [fl3|| , the present phase transition involves unphysical 
properties, including negative Z, negative specific heat in the low-temperature phase, and non-existence 
of a thermodynamic limit independent of boundary conditions in the low-temperature phase. The general 
existence of such pathologies was noted in |0J . As a increases further to 1 , B q contracts in to a point at the 
origin, q = 0. 



_2 I ■ ' ■ l_ 

1 2 

Re(q) 



Figure 9: Locus B q ({L}) for the n — > 00 limit of the L v = 2 cyclic or Mobius triangular strip for a = 0, i.e., the T = limit 
of the Potts antiferromagnet. Chromatic zeros for m = 20 (so that Z is a polynomial of degree n = 40 in q) are shown for 
comparison. 



3.3.2 Ferromagnetic range a > 1 

For the Potts ferromagnet, as T decreases from infinity, i.e. a increases above 1, the locus B q forms a lima 



bean shaped curve shown for a typical value, a = 2, in Fig. 13. Besides the generally present crossing at 



<7 = 0, the point q c ({L}) at which B q crosses the real q axis now occurs at negative q values. As was true 
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2 



_2 I ■ ' ' l_ 

1 2 



Re(q) 

Fi gure 10; Locus Bq({Ly): same as Fig. El for a = 0.1 (finite-temperature Potts antiferromagnet) . 



2 

Re(q) 



Fi gure Hi Locus Z3g({L}): same as Fig. ^ for a = 0.2. 
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Im(q) o - 



0.5 1 1.5 2 

Re(q) n 
Figure 12: Locus B q ({L}): same as Fig. |s| for a = 0.5. 

of the model on the analogous width L y — 2 cyclic and Mobius strips of the square lattice Jl3| , for physical 
temperatures, the locus B q for the Potts ferromagnet does not cross the positive real q axis. Note that this 
locus does have some support in the Re(q) > half plane, away from the real axis, which was also true of 
the analogous locus for the L y — 2 cyclic and Mobius square strip. 

3.4 B U ({L}) for Fixed q 

We next proceed to the slices of B in the plane defined by the temperature Boltzmann variable it, for given 
values of q, starting with large q. In the limit q — > oo, the locus B U ({L}) is reduced to 0. This follows 
because for large q, there is only a single dominant Xj, namely 

Al, 5 ~ q 2 + Mv + O(l) as q — > oo . (3.4.1) 

Note that in this case, one gets the same result whether one takes q — > oo first and then n = 2m — > oo, or 
n — > oo and then q — > oo, so that these limits commute as regards the determination of B u . 

We first consider values of q ^ 0, 1,2, so that no noncommutativity occurs, and (B u ) nq = {B u ) qn = B u . 
As discussed above, it is convenient to use the u plane since B u is compact in this plane, except for the cases 
q = 2 and 9 = 3, whereas B u is noncompact because of the antiferromagnetic zero-temperature critical point 
at a — l/u = 0. 

Extending the discussion in jl^] to the case of the strip of the triangular lattice, we observe that the 
property that the singular locus B u passes through the T = point u = for the Potts model with periodic 
boundary but not with free boundary conditions means that the use of periodic boundary conditions yields 
a singular locus that manifestly incorporates the zero-temperature critical point, while this is not manifest 
in B u when calculated using free boundary conditions. 
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Figure 13: Locus Bq for the n — + oo limit of the cyclic or Mobius triangular strip, {L} or {ML}, with a = 2 (finite-temperature 
Potts ferromagnet). 



For q = 10, the locus B u is shown in Fig. [L4[ Six curves on B u run into the origin, u = 0, at the angles 



given in general in eq. ( 3.2.18 ). The six corresponding complex-temperature phases contiguous to the origin, 
u = exhaust the totality of such phases; i.e., there are no such phases that are disjoint from u = 0. The 
Alj's that are dominant in these phases were given above, after eq. (3.2.18). As is evident in Fig. |l4], part 
of B u forms a approximately circular curve. The locus B u also includes a line segment on the negative real 
u axis along which Xl,5 and Xl,6 are dominant and are equal in magnitude as complex conjugates of each 
other. 

For q = 2, the locus (B u ) nq is shown in Fig. [l5| One sees that, in addition to the six curves intersecting 
at the ferromagnetic zero-temperature critical point u = 0, B u has multiple points at u = ±z where four 
curves intersect. The complex-temperature phases in the vicinity of u = were determined above after eq. 
( |3 .2.18 ), and these exhaust all complex-temperature phases, i.e. there is none that does not extend in to 
u = 0. 

For g = 3, the locus B u was given as Fig. 10 in Ref. [[l2| and is shown with associated partition 
function zeros in Fig. |l7]. In this case, in addition to the six phases that are contiguous at u = 0, there is 
evidently another O phase that includes the negative real axis for u < — 1, and a complex conjugate pair of 
O phases extending outward from the intersection points at u = e ±27 ™/ 3 toward the upper and lower left. 
At these intersection (multiple) points, six curves on B u intersect, just as was true at u = 0. The existence 
of the intersection points on B u at u = e ±27 ™/ 3 for the strip would suggest that such points could also be 
present on the locus B u for the q = 3 Potts model on the full triangular lattice. Since this model has not 
been exactly solved, one can only try to infer B u from complex-temperature (Fisher) zeros of the partition 
function calculated for large triangular lattices 45, |o|, EH]]. These are consistent with this possibility (see, 
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e.g., Fig. 1 of 0] or Figs. 5-7 of JsT]) although the zeros show such a high degree of scatter in the Re(a) < 
half-plane that one cannot draw a very decisive conclusion from them. 
For q = 4 we show the locus B u in Fig. [l^. In this case, we remark, in 

For the study of the zero-temperature critical point of the Potts antiferromagnet on the L y — 2 cyclic 
and Mobius strips of the triangular lattice, it is useful to display the singular locus B a in the a plane, since 
the critical point occurs at a = 0. We show these plots for q = 2 and 3 in Figs. |l6| and |l||. For q = 2, the 
boundary B a has a multiple point at a — ±i where four branches intersect. For q = 3, four curves on B a 
pass through a = at the angles given in eqs. ( 3.2.28 ) and ( 3.2.29 ). Note that that the curves that leave 
the origin at the angles 4> given in eq. ( 3.2.28| ) rapidly bend toward the vertical and then back toward the 
left, so that the complex-temperature phase that is contiguous with u = and lies in the angular wedges 
between 76° and 150°, and its complex-conjugate, are rather narrow. The multiple points at a = e ±27 ™/3 
have been described above. 



In Figs. |21| and |20| we show the analogous loci B u for the L y = 2 cyclic or Mobius strip of the square 
lattice studied in |l3| for comparison. A particularly striking difference is that, since the q = 3 Potts 
antiferromagnet is (is not) critical on the L y = 2 strip of the triangular (square) lattice, the resultant locus 
B u passes through (does not pass through) 1/u = a = 0, respectively. The Ising model, q = 2 has both 
ferromagnetic and antiferromagnetic T — critical points on both the L y = 2 square and triangular lattice 
strips, and hence for both strips, B u passes through 1/u = a = 0. 
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Figure 14: Locus B U ({L}) for the n — > oo limit of the cyclic or Mobius triangular strip, {L} or {ML}, with q = 10. Partition 
function zeros are shown for m = 20, so that Z is a polynomial of degree e = Am = 80 in v and hence, up to an overall factor, 
in u). 



3.5 Connections Between B for Strips and 2D Lattices 

In earlier work [p^ , |54l fi"3"| , it was shown that although the physical thermodynamic properties of a discrete 
spin model are, in general, different in ID or infinite-length, finite-width strips, which arc quasi-lD systems, 
and in higher dimensions, nevertheless, exact solutions for B u in ID and quasi- ID systems can give insight 
into B u in 2D. This was shown, in particular, for the q = 2 Ising special case of the Potts model, where the 
comparison can be made rigorously since the model is exactly solvable in 2D. It will often be convenient 
below to use the equivalent locus in the a — 1/u plane, B a . In |2q, ^5|, [l2| Q it was also noted how exact 
calculations of B q and W on infinite-length, finite-width strips can give information about the behavior of 
these quantities on 2D lattices. 
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Figure 15: Locus B U ({L}): same as Fig. ^ for q = 2. 
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Fi gure 16i Locus B a for the n — ► oo limit of the cyclic or Mobius triangular strip, {-£/} or {ML}, with q = 2. 
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Figure 17: Locus B U ({L}): same as Fig. ^ for q = 3. 
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Fi ffure 18: Same as Fig. Ilq for q = 3. 
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Fi gUTC 20i Locus 3 U — $a for the L y = 2 square strip with q = 2. Partition function zeros are shown for m = 20 (e = 60). 
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Figure 21: Locus B u for the L y = 2 square strip with q = 3. Partition function zeros are shown for m = 20 (so that Z is a 
polynomial of degree e = 3m. = 60 in v and hence, up to an overall factor, in u). 
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Fi glire 22: Locus 3 a for the L y = 2 square strip with g = 3. Partition function zeros arc shown for m = 20 (so that Z is a 
polynomial of degree e = 3m = 60 in v and hence a). 
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In Ref. jy|, the connection between B a on strips and in on the full 2D square lattice was studied; here 
we shall do this for the triangular lattice. Again, it is natural to start with the q = 2 Ising case, where 
one has exact results both for the strip and the full triangular lattice. The comparisons made in |T^ | and 
here enable one to formulate a reasonably systematic procedure for making transformations on an exactly 
calculated locus B a for the Potts model on an infinite-length, finite-width strip of a given lattice, in order to 
construct at least a qualitatively correct locus B a on the corresponding 2D lattice. As emphasized in [ O , 
this represents a new and powerful way of obtaining information about complex-temperature phase diagrams 
of spin models in 2D (indeed perhaps also higher dimensions) from exact results on strips. It will be recalled 
that the conventional approach has been the rather laborious procedure of performing exact calculations 
of the partition function on finite 2D lattices of various sizes with various boundary conditions p5[-[p9|. 
Although some features of B a for the Potts model, such as the circle \a— 1| = ^/q for Re(a) > are evident 
with this conventional procedure, the comlex-temperature (Fisher) zeros often exhibit such considerable 
scatter for Re (a) < that it is difficult to infer the structure of the locus B a in this region. One method was 
to combine calculations of the zeros with dlog Pade and differential approximant analyses of low-temperature 
series expansions so as to localize accurately certain points on the complex-temperature boundary Jf9[-|3l|. 
This method was motivated by the fact that complex-temperature singularities in such quantities as specific 
heat and magnetization for the 2D Ising model can be calculated exactly and they occur at known points on 
the complex-temperature boundaries for respective 2D lattices; it was also shown that this was true for the 



susceptibility, where no exact calculation exists 139, 52 . However, while this method avoids the problem of 



the scatter in Fisher zeros, it only localizes certain points on B a . The present method is complementary in 
that it enables one to gain, at least qualitatively, an idea of the global structure of B a . 

We proceed with our exact comparison for the q = 2 Ising case and start by recalling the situation for the 
square lattice The locus B u for the square lattice consists of the well-known union of circles \u — 1| = y/2 
and \u + 1| = \/2-2 The procedure for transforming the locus B u found for the infinite-length, width L y = 2 
strip in Jl3| (see Fig. 20) involves two main steps. First, one retracts each of the curves going through the 
origin u = so that they no longer pass through this point. It is clear that this retracting is necessary not 
just for q — 2 but for the general g-state Potts ferromagnet since the (reduced) free energy and magnetization 
of this model are analytic in the neighborhood of u = 0, i.e., they have low-temperature series expansions 
a finite radius of convergence. Given the inversion symmetry B u = B a that holds for the square lattice and 
strips of this lattice (but not for the triangular lattice), this retracting means that the curves are also pulled 
back from the origin of the a plane. The second step in the procedure is to incorporate the feature that the 
2D Potts ferromagnet has a finite-temperature phase transition. To build in this feature, one takes the two 
complex-conjugate ends of the curves with Re(u) small and positive that have been pulled back from the 
origin u — and connects them so that they cross the positive real u axis at a point upm-fm in the interval 
< u < 1; by the inversion symmetry this has the effect of connecting the two other ends in the Re(u) > 
half plane and having them cross the real axis at upm-afm — u pm-fa/- The third step is to build in the 
feature (3^, ^2j that if a lattice has an even coordination number, then the Ising model boundary B u is 
symmetric under u — > — it; in the present case, the locus B u for the L y — 2 strip does not have this property 
because the coordination number of the vertices on the strip is 3, while the locus B u for the square lattice 
does have the property. We thus connect the curves in the left-hand half plane Re(u) < in such a way 

3 In Refs. fcisj, |52|] the variable z = e~ 2jfjsi "g = e~ Kp ° tts was used for what we denote as u here and the variable u was 
used for e~ 4Klsin s = F °" s , i.e. for what would be denoted v? here. Thus, the circles \u ± 1| = \/2 map to the limagon of 

Pascal in the u 2 plane [B9[ . 
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as to have this property. As discussed in 13 , the intersection points at u — ±i are the same for both the 
strip and the square lattice, and, indeed, on wider strips. Evidently, each step of this procedure is based 
on fundamental principles, and, as long as one limits oneself to obtaining the qualitative features of the 
locus B u for the 2D case from the L y = 2 strip, nothing is ad hoc. Of course, one cannot predict the actual 
values of the critical point Upm-fm = u p^m-afM' but ^ ms * s was n °t the § oa l; there are powerful ways of 
determining the critical point via series analyses even in cases where a spin model cannot be solved exactly in 
2D. (Note that for the purpose of obtaining information about B u for a model in 2D, it can be advantageous 
to use strips with width L y = 2 rather than wider strips, since the locus B u becomes progressively more 
complicated as L y increases, with more curves passing through u = p3[.) 

Next, we demonstrate the corresponding comparison with the locus B u for the q — 2 Ising special case 
of our new results for the Potts model on the L y = 2 strip of the triangular lattice and B u for the model on 
the full 2D triangular lattice. Our calculations are shown in Figs. |l5| and [l6|. For the Ising model on the full 
2D triangular lattice, B u consists of the union of an oval curve that crosses the real (imaginary) u axis at 
±z with the line segment 1/V3 < Im{u) < oo and its complex conjugate on the imaginary u axis [5^ ]. Note 
that, just as was true in the case of the square strips and square lattice (and the honeycomb lattice [p2|), 
the locus B u for the Ising model on both the current strip and the full 2D triangular lattice has intersection 
points at u = ±z. (In terms of the variable denoted u 2 in our current notation (and u in the notation of 
|j52f), the Ising locus for the triangular lattice is the union of the semi- infinite line segment — oo < u 2 < — 1 
and the circle \u 2 + (1/3) | = 2/3, and the above two intersection points map to the single intersection point 
at u 2 — —1.) Following the same procedure as for the square strip, the first step is to pull back each of the 
curves (of which there are now six rather than four) from the origin. On the imaginary axis, this retracting 
produces two complex-conjugate semi-infinite line segments, while for the other four curves, it allows one 
to connect them smoothly, via step 2, to form the oval. Step 3 is not necessary here since both the cyclic 
L y = 2 strip of the triangular lattice and the full triangular lattice have the property that each vertex has 
even coordination number (A = 4 for the strip and A = 6 for the 2D lattice) so that B u is invariant under 
the replacement u — > — u |32| |55| . 

It is also interesting to observe that the similarities between the Ising model on the current strip and on 
the full 2D lattice are much stronger for the antiferromagnet than for the ferromagnet. In contrast to the 
case of the Ising ferromagnet, which has a finite-temperature critical point on 2D lattices but only a zero- 
temperature critical point on ID and quasi-lD lattices, the Ising antiferromagnet has a zero-temperature 
critical point on both the full 2D triangular lattice [Q and on the infinite-length, width L y — 2 strip of the 
triangular lattice, as a consequence of the frustration in this model. Thus in this case, there are not just 
similarities in the complex-temperature properties of the model, but in the actual physical thermodynamics 
also. Concerning the complex-temperature phase diagram, in the vicinity of this zero-temperature critical 
point, for the full 2D triangular lattice, B a is just the inverse image of the locus B u given above, namely the 
union of an oval intersecting the real (imaginary) a axes at ±-\/3 and ±i, and a finite line segment on the 
imaginary axis stretching from a — — v3i to a = \/3i. In the neighborhood of the origin, a = 0, the locus 
B a for the present strip is exactly the same, as can be seen in Fig. |l^. 

Having shown the correspondences between the complex-temperature boundaries B u and the associated 
complex-temperature phase diagrams for the exactly solved Ising case, we now use this result as a tool to 
suggest features of B for the 2D Potts model for other values of q, where it has not been exactly solved. 
Since the complex-temperature zeros have usually been presented in the a plane, we follow this convention 
here. The case q — 3 on strips of the square lattice and on the full square lattice was discussed in jl3|. In 
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Figs. |2lj and |22| wc show B u and B a for the L y — 2 strip of the square lattice. First, the fact that B a crosses 
the real axis at a = —2 for the strip is the same as is true for the square lattice jlj|; in the latter case, this 
follows from the existence of the zero-temperature critical point for the q — 3 antiferromagnet, the resultant 
fact that the singular locus B a passes through a — 0, and the duality of the model, according to which if a 
point a is on B a , then the dual image a c [ = (q + a — I) /(a — 1) is also on B a . Second, the existence of the 
intersection points at a = e ±2 ™/ 3 on B a for the strip suggests that these also occur on B a for the square 
lattice. We now extend these two comments to explore the global structure of B a . Since we know that the 
2D q = 3 Potts model on the square lattice is analytic in the neighborhood of the T = ferromagnetic point 
u = 0, our first step is to pull back the four curves that pass through u = in Fig. ^l] (equivalently, off to 
infinity in Fig. p2| ). Since we know that the 2D model has a finite-temperature critical point (the exact value, 
a = l|s<7ri3 is not crucial here), we connect the two ends of the curves in the Re{a) > half-plane together 
in Fig. ^2] so that they cross the positive real a plane, making the ferromagnetic critical point clpm-fm > 1- 
The ends of the two curves pointing in the upper left and lower left directions in the Re(a) < half-plane 
are left as is, in accordance with the finding from series analyses |l9| that there are such singularities at 
a e = —1.71(1) ± lA6(q)i, consistent with lying on endpoints of curves on B a . Since we know that the q = 3 
Potts antiferromagnet has a T — critical point on the square lattice, and hence B a passes through a = 0, 
we move the four curves that pass through the point a = —1/2 (forming an intersection point) over to the 
right, so that at least one branch passes through a = 0; the exact solution for the strip suggests that B a 
for the square lattice may also have an intersection (multiple) point at a = Q. The complex-temperature 
zeros for the q = 3 Potts model that have been calculated for patches of the square lattice ^||, |Z(], ^] 
are consistent with this suggestion (see, e.g., Fig. 1 of j|9|]). A third suggestion is that the intersection 
points on B a at a ~ —1.85 ± 0.85z for the strip have analogues on the locus B a for the full square lattice. 
This suggestion is consistent with the patterns of complex-temperature zeros that have been calculated for 
patches of the square lattice, but the scatter is too great to provide a strong test. 

We next use our new results on the q — 3 Potts model on the L y — 2 strip of the triangular lattice. 

3.6 Thermodynamics of the Potts Model on the L y = 2 Strip of the Triangular 
Lattice 

3.6.1 Ferromagnetic Case 



The Potts ferromagnet (with real q > 0) on an arbitrary graph has v > so, as is clear from eq. (1.5), 
the partition function satisfies the constraint of positivity. In contrast, the specific heat C is positive for 
the model on the (infinite-length limit of the) L y = 2 triangular strip is positive if and only if q > 1 (for 
any choice of longitudinal boundary conditions). For q = 1, f nq = 2 In a = 2K and C vanishes identically. 
Since a negative specific heat is unphysical, we therefore restrict to real q > 1. For general q in this range, 



the reduced free energy is given for all temperatures by / = (l/2)lnAs,i as in (2.12) (independent of the 
different longitudinal boundary conditions, as must be true for the thermodynamic limit to exist). Recall 
that As,i = It is straightforward to obtain the internal energy U and specific heat from this free energy; 

since the expressions are somewhat complicated, we do not list them here. We show a plot of the specific 
heat (with fc^ = 1) in Fig. ^3|. One can observe that the value of the maximum is a monotonically increasing 
function of q. 
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Figure 23: Specific heat (with ks = 1) for the Potts ferromagnet on the infinite- length, width L y = 2 strip of the triangular 
lattice, as a function of K = J/(kgT). Going from bottom to top in order of the heights of the maxima, the curves are for 
9 = 2,3,4,6. 



The high-temperature expansion of U is 



q 



i+(±—Hk + o(k 2 ) 



For the specific heat we have 
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The low-temperature expansions (K — > oo) arc 
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and 
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as K -> oo (3.6.4) 



In general, the ratio p of the largest subdominant to the dominant Xj 's determines the asymptotic decay 
of the connected spin-spin correlation function and hence the correlation length 

1 



In p 



Since Al,5 and Al,2 are the dominant and leading subdominant A/s, respectively, we have 

Ai,2 



Pfm 



X 



L,5 



(3.6.5) 



(3.6.6) 



and hence for the ferromagnetic zero-temperature critical point we find that the correlation length diverges, 
as T — > 0, as 

£,fm - (2g) _1 e 3K , as K -> oo (3.6.7) 
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3.6.2 Antiferromagnetic Case 

In this section we first restrict to the real range q > 3 and the additional integer value q = 2 (Ising case) 
where the Potts antifcrromagnet exhibits physically acceptable behavior and then consider the remaining 
interval < q < 3 where (except for the trivial f nq for q = 1) it exhibits unphysical properties. For q > 3, 



the free energy is given for all temperatures by ( 2.12 ), as in the ferromagnetic case but with J negative, and 
is the same independent of the different longitudinal boundary conditions, as is necessary for there to exist 
a thermodynamic limit. 

We show plots of the specific heat, for several values of q, for the Potts antiferromagnet on the (infinite- 
length limit of the) L y — 2 strip of the triangular lattice in Fig ^4|. In contrast to the ferromagnetic case, 
the maxima of C do not increase monotonically with q; they occur at 0.32 for q = 2 and 0.79 for q = 3, 
after which the values of the maxima decrease with increasing integral q (e.g., they are 0.45 and 0.24 for 
q = 4 and q = 6). The fact that the curve for C for the q = 2 (Ising) antifcrromagnet exhibits a very broad 
maximum can be inferred to be a consequence of the frustration that is present in this model. 




Figure 24: Specific heat (with fcg = 1) for the Potts antiferromagnet on the infinite-length, width L y = 2 strip of the 
triangular lattice, as a function of —K = —J/(kgT). Going downward in order of the heights of the maxima, the curves are 
for q = 3,4,2,6. 



The high-temperature expansions of U and C are given by ( |3.6.1| ) and ( |3.6.2| ); more generally, these 
expansions also apply in the range < q < 3. As discussed above, the Ising case q = 2 is one of the 
cases where one must take account of noncommutativity in the definition of the free energy and hence of 
thermodynamic quantities. If one sets q — 2 first and then takes n — > oo, then / = f nq — (1/2) In Al.5| 9 =2 



where \L,a\q — 2 was given in eq. (3.2.7), and the low-temperature expansions are 
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(3.6.8) 



(3.6.9) 
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For the range q > 3, the low-temperature expansions are given by 
{-J)e K 



and 



U 



C = 
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(3.6.10) 



(3.6.11) 



{q-2f 

Note that for the antiferromagnetic case, U(T = 0) = for q > 3, but U(T = 0) = - J/2 = +\J\/2 for g = 2. 
The vanishing value of U at T = for q > 3 means that the Potts model can achieve its preferred ground 
state for this range of q, while the nonzero value of U(T = 0) for the Ising antiferromagnet is a consequence 
of the frustration that is present in this case. Similarly, the fact that the specific heat vanishes less rapidly 
for the Ising antiferromagnet than for the Potts model with q > 3 reflects the frustration that is present 
in the Ising case, which is not present for q > 3. Note that the apparent divergences that occur as q — > 2 
m eqs. (E.6.1CQ and ([3. 6. IIP are not actually reached here since these expressions apply only in the region 
q > 3 (the discrete integral case q = 2 was dealt with above). 

For the zero-temperature critical points in the q = 2 and q = 3 Potts antiferromagnet, 



PAFM,q=2, 3 



A 



L,5 



(3.6.12) 



Using the respective expansions ( |3.2.19| )-( |3.2.20] ) and ( |3.2.21 )-( 3.2.24 ), we find that the correlation lengths 
defined as in (3.6.5) diverges, as T — > 0, as 



AFM,q=2 



and 



-K/2 



^AFM, q =3 ~ e K , as if -> -oo (3.6.14) 

Next, we consider the range of real nonintegral < q < 3. The first pathology is that the Potts 
antiferromagnet on the infinite-length limit of the L y = 2 triangular strip, defined with case of cyclic or 



as K — ► — oo 



(3.6.13) 



Mobius longitudinal boundary conditions has a phase transition at the temperature T p .l given in eq. ( 3.3.7 ) 
for < q < 3 (except for f nq for the integral valuesg — 1,2), while, in contrast, if one uses free boundary 
conditions, then there is no phase transition at this temperature and, although there is a phase transition 
for < q < 2, it occurs at the temperature T p $ given in eq. (2.2.6), which, in general, is not equal to 
T P: l. It follows that there is no well-defined thermodynamic limit for the Potts model with < q < 3 
and q ^ 1,2. The Ising case q = 2 has been dealt with in the preceding subsection. Concerning the value 
q = 1, as discussed earlier, one encounters noncommutativity in defining the free energy. If one takes q = 1 
to start with and then n — > oo, the thermodynamic limit does exist, independent of boundary conditions, 
and / = f nq — 2K, U = — 2J = 2\J\, and C = 0. If one starts with q ^ 1, takes n — * oo, calculates 
f qn , and then takes q — > 1, the thermodynamic limit does not exist since the result differs depending 
on whether one uses free longitudinal boundary conditions or cyclic (equivalently Mobius) longitudinal 
boundary conditions. Specifically, for this single value q = 1, T Pi s is equal to T P: l, having the value given by 
k B T p , s = J/ln[(l/2)(-l + s/5)] ~ 2.078| J|; in the high-temperature phase, T > T P:S , f qn = (1/2) In X L>5 , 
independent of longitudinal boundary conditions, but in the low-temperature phase, T < T Pj s, the expression 
for fqn is different for the open and cyclic (equivalently, Mobius) strips. There are also other unphysical 
properties, such as a negative specific heat and a negative partition function for certain ranges of temperature. 
These is similar to what was found for the analogous strip of the square lattice. 
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4 Summary 



In this paper we have presented exact calculations of the partition function Z of the q-state Potts model 
and its generalization to real q, the Potts model, for arbitrary temperature on n-vertex strip graphs, of 
width L y = 2, of the triangular lattice with free, cyclic, and Mobius longitudinal boundary conditions. 
These partition functions are equivalent to Tutte/Whitney polynomials for these graphs. The free energy is 
calculated exactly for the infinite-length limit of these ladder graphs and the thermodynamics is discussed. 
Considering the full generalization to arbitrary complex q and temperature, we determine the singular locus 
B in the corresponding C 2 space, arising as the accumulation set of partition function zeros as n — > oo. 
In particular, we study the connection with the T = limit of the Potts antiferromagnet where B reduces 
to the accumulation set of chromatic zeros. Comparisons are made with our previous exact calculation of 
Potts model partition functions for the corresponding strips of the square lattice. Our present calculations 
yield, as special cases, several quantities of graph-theoretic interest, such as the number of spanning trees, 
spanning forests, etc., which we record. 

Acknowledgment: The research of R. S. was supported in part at Stony Brook by the U. S. NSF grant 
PHY-97-22101 and at Brookhaven by the U.S. DOE contract DE-AC02-98CH10886.(] 

5 Appendix 
5.1 General 

The formulas relating the Potts model partition function Z(G, q, v) and the Tutte polynomial T(G, x, y) were 
given in |l3) and hence we shall be brief here. The Tutte polynomial of G, T(G, x, y), is given by |6)-§ 

T(G,x,y)= J2 (x-l) k ^- k ^(y-lY lG,) (5.1.1) 

G'CG 

where k(G'), e(G'), and n(G') = n(G) denote the number of components, edges, and vertices of G', and 

c(G') = e(G') + fc(G') - n(G') (5.1.2) 

is the number of independent circuits in G'. As stated in the text, k(G) = 1 for the graphs of interest here. 
Now let 

x = l + - (5.1.3) 

v 

and 

y = a = v + l (5.1.4) 

so that 

q=(x-l)(y-l) . (5.1.5) 

Then 

Z(G,q,v) = (x-l) k ^(y-ir^T(G,x,y) . (5.1.6) 
For a planar graph G the Tutte polynomial satisfies the duality relation 

T(G,x,y)=T(G*,y,x) (5.1.7) 



4 Accordingly, the U.S. government retains a non-exclusive royalty-free license to publish or reproduce the published form of 
this contribution or to allow others to do so for U.S. government purposes. 
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where G* is the (planar) dual to G. As discussed in |13[ , the Tutte polynomial for recursively defined graphs 
comprised of m repetitions of some subgraph has the form 



T(G m ,x,y) = y^Qr ]G j(AT,G,j) T 



(5.1.8) 



There are several special cases of the Tutte polynomial that are of interest. One that we have analyzed 

in the text and in previous papers is the chromatic polynomial P(G,q). This is obtained by setting y = 0, 

i.e., v = — 1, so that x = 1 — q; the correspondence is P(G,q) = (—q) k ^(—l) n T(G,l — q,0). A second 

special case is the flow polynomial J9|, [II], |llj F(G, q), obtained by setting x = and y = 1 — q: F(G, q) = 
(_l)e(G)-n(G)+fe(G) T ( G)0)1 _ q j For planar Gj 

given the relation ( J5 . 1 . 7| ) , the flow polynomial is, up to a 
power of q, proportional to the chromatic polynomial: F{G,q) oc P(G*,q). A third special case is the 
reliability polynomial fTo|| , which will not be considered here. 

5.2 Triangular Strip with Free Longitudinal Boundary Conditions 

The generating function representation for the Tutte polynomial for the open strip of the triangular lattice 
comprised of m + 1 squares with edges joining the lower-left to upper right vertices of each square, denoted 
S m , is 



T T (S m ,x,y;z) = T(S m ,x,y)i 



m=0 



We have 

where 
with 

and 



r T (S,x,y;z) 



Af T (S,x,y;z) 



V T {S,x,y;z) 
N T {S, x, y; z) = A TiS ,o + A T . s ,iz 

A T ,s,o = x(x + l) 2 + 2xy + y(y + 1) 
A T ,s,i = ~x 3 y 2 



(5.2.1) 

(5.2.2) 

(5.2.3) 

(5.2.4) 
(5.2.5) 



T> T (S,x,y,z) = 



1- [{x + lf +y{y + 2)]z + x 2 y 2 z 2 

2 

H(1-Xt,sjz) 

3=1 



with 



where 



A 



T,S,(1,2) 



(a: + I) 2 +y{y + 2)±(l + x + y)y/ R T 



R T = ( X + \f + (y + I) 2 - 1 - 2xy . 
The corresponding closed-form expression is given by the general formula |1J 



T(S m ,x,y) 



A 



T.S.1 — At,S.2 



(At,s,i)'™ + 



At,S,0^T,S,2 + At,s,i 



A 



T,S,2 



A 



T.S.I 



(AT,S,2) r 



(5.2.6) 

(5.2.7) 
(5.2.8) 

(5.2.9) 



It is easily checked that this is a symmetric function of the As,,-, j = 1, 2. 
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5.3 Cyclic and Mobius Strips of the Triangular Lattice 

We have 

6 

T(L m ,x,y) =J2 c T,LA*T,L,j) m (5.3.1) 
3=1 

and 

6 

T(ML m , x,y) = Y^ c T ,ML,j(\TML, 3 ) m (5.3.2) 

3=1 

It is convenient to extract a common factor from the coefficients: 

c T Gi= < ^ ± , G = L,ML . (5.3.3) 
x — 1 

Of course, although the individual terms contributing to the Tutte polynomial are thus rational functions of 

x rather than polynomials in x, the full Tutte polynomial is a polynomial in both x and y. We have 

Ar.Mij = A t ,lj , 3 = 1, -,6 (5.3.4) 

Xt,l,i = 1 (5.3.5) 

and 

At,5 = At,s,i j At,6 = At,s,2 (5.3.6) 

where \t,s,i an d At,s,2 were given in eq. ( 5.2.7 ) with fl5.2.gf) . The Xr,L,j, j = 2,3,4, are solutions of the 
cubic equation 

e - (3 + 2x + 2y + y 2 )e + [(a: + l) 2 + 2xy(y + 1)]^ - x 2 y 2 - (5.3.7) 

We label these At,Lj's in a manner corresponding to the A^j's occurring in the Potts model partition 
function, so that At,l,5 is dominant in the region in (x, y) equivalent to the PM phase, and so forth for the 
others. We note that 

^t,l,2^t,l,3^t,la = At, l, 5 At, l, 6 = x 2 y 2 . (5.3.8) 
The coefficients for the cyclic strip are 

Qt,l,i = [(x - l)(y - I)} 2 - 3(x - l)(y - 1) + 1 (5.3.9) 

cr,L,j = xy - x — y for j = 2, 3,4 (5.3.10) 

and 

ct,lj = 1 for j = 5,6. (5.3.11) 

These are symmetric under interchange of x <-> y, which is a consequence of the fact that the clj are 
functions only of q, which, by eq. (^.1.5 ) is a symmetric function of x and y. For the Mobius strip, we find 

c t ,ml,i = -1 (5.3.12) 

and 

CT,ML,5 = C T ,MLfi = 1 • (5.3.13) 
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The ct,mlj, j — 2,3,4 can be calculated from the generating function using the formula given in [p~3| - 
Their form is more complicated, as is evident from the results given in |^2| for the corresponding chromatic 
polynomial for this case. We write the generating function as 

oo 

Ft (ML, x, y,z)=J2 T(ML m , x, y)z m - 1 . (5.3.14) 

m— 1 

with 

T T (ML, x, y, z) = (5.3.15) 

V T (ML,x,y, z) 

6 

V T {ML, x, y, z) = V T (L, x, y, z) = J](l - X L>j z) (5.3.16) 

5 

M T (ML, x, y,z) = Y^ A ML<j z j (5.3.17) 
j=o 

A M Lfi = v(x + y + y 2 ) (5.3.18) 
A M l,i = yx- Ay 2 x - Ay 3 - 3y 4 - 5y 3 x - y 5 - 2x 2 y - x 3 y - x 2 y 2 - x 2 y 3 + 2y + 2x + 3x 2 + x 3 (5.3.19) 

A M l,2 = y{-2x 2 y + Ayx + 3x + 3y + x 2 y 3 -2x 3 y + &x 2 y 2 + 2xy 4 

+ Ay 3 x + y 4 x 2 + 3x 3 y 2 + by 2 + 3y 3 + y 4 + 6y 2 x + 2x 2 ) (5.3.20) 



where 



and 



with 



A M l,3 = -y{-2x 2 y - 3yx - x - y - x 4 y + 3x 2 y 3 - 2x 3 y + 8x 2 y 2 + 2xy 4 



+ 2x 3 y 4 - x 3 y 3 + 5y 3 x - x 3 + x 4 y 2 + 2y 4 x 2 + Ax 3 y 2 - 2y 2 



y 3 + 3y z x - 2x z ) 



(5.3.21) 



Amla = x A y\-l + x V - x V -x z +y z -2x + y + 2y A x + x l y + y 3 + yx) 



(5.3.22) 



and 



A M l,5 = ~y 4 x 4 (y - 1) . 



(5.3.23) 
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5.4 Special Values of Tutte Polynomials for Strips of the Triangular Lattice 

For a given graph G — (V,E), at certain special values of the arguments x and y, the Tutte polynomial 
T(G,x,y) yields quantities of basic graph-theoretic interest (|]-@- We recall some definitions: a spanning 
subgraph was defined at the beginning of the paper; a tree is a connected graph with no cycles; a forest is 
a graph containing one or more trees; and a spanning tree is a spanning subgraph that is a tree. We recall 
that the graphs G that we consider are connected. Then the number of spanning trees of G, Nst(G), is 



N ST (G)=T(G,l,l) , 
the number of spanning forests of G, Nsf(G), is 

N SF (G)=T(G,2,l) , 
the number of connected spanning subgraphs of G, Ncssg(G), is 

N C ssg(G) = T{G,1,2) , 
and the number of spanning subgraphs of G, Nssg(G), is 

N S sg(G) =T(G,2,2) . 



From the duality relation (5.1.7), one has, for planar graphs G and their planar duals D* , 

N ST (G) = N ST (G*) , N SSG (G) = N SSG (G*) 

and 

N SF (G) = N C ssg{G*) . 
From our calculations of Tutte polynomials, we find that 



9^5 \ (7 - 3V5 



(5.4.1) 

(5.4.2) 

(5.4.3) 

(5.4.4) 

(5.4.5) 
(5.4.6) 

(5.4.7) 



17\/2 \ r — " <» 

N SF (S m ) = Ncssam = (l2 + J [2(3 + 2^2 )J + (l2 



17^2 



2 ) 



2(3 - 2^2) 



N C ss G {S m ) = N SF (S*J = (7 + 5V2 ) [2(3 + 2V2 )1 + (7 - 5V2 ) [2(3 - 2V2 ) 



NssG(S m ) = N SSG (S* m ) = 2 4m + 5 



(5.4.8) 



(5.4.9) 



(5.4.10) 



For the cyclic L y = 2 strip of the triangular lattice, L m , we first note that for m > 3, L m is a (proper) 
graph, but for m = 1 and m = 2, L m is not a proper graph, but instead, is a multigraph, with multiple 
edges (and, for m — 1, loops). The following formulas apply for all m > 1: 



2?n 

N ST (L m ) = N ST (L* m ) = — 



7 + 3V5 



7-3V5 



(5.4.11) 
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r -i rri r 

N SF (L m ) = N C ssa(L*J = 1 - ^ (A Tj ) m + [2(3 + 2^2 )J +[2(3-2^2) 

3=2,3,4 



where Xt,L,j, j = 2,3,4 are the roots of eq. (5.3.7) for x = 2,y = 1, viz., £ d - 10£ 2 + 17£ - 4 = 0: 



>^T,L,2(x 


= 2,y = 


1) 


= 0.280176 


^T,L,3(x 


= 2,y = 


1) 


= 1.803442 


Xt,l,a(x 


= 2,y = 


1) 


= 7.916382 



(5.4.12) 

(5.4.13) 

(5.4.14) 
(5.4.15) 



We also calculate 

Ncss G (L m ) = N SF (L* r , 



6m, 
"(2+—) 



1 + m 



1 + m 
(3 + V2 ) 



(3-^2) 



2(3 + 2V2 



7 



2(3 - 2V2 ) 



(5.4.16) 



N SSG (L m ) = N SSG (L* m ) = 2 4m . (5.4.17) 

Since T(G m ,x,y) grows exponentially as m — > 00 for the families G m = S m and L m for (x,y) = (1, 1), 
(2,1), (1,2), and (2,2), one defines the corresponding constants 



z se t({G})= lim n(G)" 1 lniV sef (G) , set = ST, SF, CSSG, SSG 

n(G)—» co 



(5.4.18) 

where, as above, the symbol {G} denotes the limit of the graph family G as n(G) — > 00 (and the z here 



should not be confused with the auxiliary expansion variable in the generating function (5.2.1) or the Potts 
partition function Z(G, q, v).) General inequalities for these were given in ]l3|]. Our results yield 



zst({G}) = - In I 



0.962424 for G = S,L, ML 



(5.4.19) 



z SF {{G}) = z CS sg{{G}) = hn[2{3 + 2V2)}~ 1.22795 for G=S,L,ML (5.4.20) 



and 



z S sg({G}) =2 ln2~ 1.38629 for G = S,L,ML (5.4.21) 

The cyclic L y = 2 strip of the triangular lattice is a K-regular graph with /c = 4. It is therefore of interest 
to compare the value of z$t that we have obtained for the n — > 00 limit of this family with an upper bound 
for ft-regular graphs with k > 3 [p7|, pq| , viz. 



(«-!)" 



[k(/c-2)]("/ 2 )- 



2st({G}) < z S T,u.bX{G}) , where ({G}) = In 

For this purpose we define r as the ratio of the left- to the right-hand side of eq. (5.4.22). We have 

Uj I+3V1 ] 



(5.4.22) 



r ST ({L}) 



3 In 



0.791 



(5.4.23) 



5 A re-regular graph is a graph in which all of the vertices have the same degree, re. 



39 



Another comparison of interest is the ratio of zst for these L y — 2 strips with zst for the full 2D triangular 
lattice, which has the value |5(J 

z ST tri = 2z hc = — (1 - 5~ 2 + 7~ 2 - ll- 2 + 13~ 2 -...) = 1.615329736... (5.4.24) 

^#^^ 0.596 (5.4.25) 

z ST {tri) K ' 

Finally, we may also compare the z values above with those for the infinite-length limit of the L y = 2 strip 
of the square lattice (with free transverse boundary conditions). This is J5?| [lj| 

z ST = - ln(2 + V3 ) ~ 0.658479 (5.4.26) 

Further, the calculation of the Tutte polynomial for the open and/or cyclic L y = 2 strip of the square lattice 
11 yields 



z SF = - ln[2(2 + V3 )] ~ 1.00505 (5.4.27) 
I ln( 5 + ^ ) ~ 0.758832 (5.4.28) 



2 

ZCSSG = 



2 

zssg = ^ In 2= 1.03972 (5.4.29) 

Thus, one observes that each of these four quantities is greater for the (infinite-length limit of the) L y — 2 
strip of the triangular lattice than for the (same limit of the) L y = 2 strip of the square lattice. For zst, this 
follows from the fact that one can obtain the strip of the triangular lattice from that for the square lattice 
by uniformly adding a diagonal bond to each of the squares |^l| . The inequalities for the other quantities 
Z SF , zcssg, and zssg are presumably a consequence of this fact also. 

5.5 Tutte Polynomials for Dual Graphs 



Since the Tutte polynomial satisfies the duality relation ( |5 .1.7 ) for a planar graph G, our calculations of the 



Tutte polynomials T(G m ,x,y) for the open and cyclic L y = 2 strips of the triangular lattice, S m and L m , 
also yield the corresponding Tutte polynomials for the duals of these graphs. For m > 3, the dual of L m is 
a graph with n = 2m + 2 vertices comprised of the circuit graph Cim together with two additional vertices, 
which we denote v e and v such that, if we label the vertices on Cim as Vj, j = 1, .., 2m, then each of the 
Vj with j even is connected by an edges to v e and each of the Vj with j odd is connected to v a . We shall 
denote this graph as an "alternating wheel" , A2 m +2 ■ Then 

T(A 2m+2 , x, y) = T(L m , y, x) , for m>3 (5.5.1) 

The cases m — 1 and m — 2 can be dealt with in a similar manner. 

The dual of the open strip graph S m is the multigraph formed by a line of 2(m + 1) vertices, of which 
the interior 2m vertices are connected by single edges to a single additional vertex, and the two end vertices 
on the line are connected to this additional vertex by double edges. We have T((S m )* ,x,y) — T(S m ,y,x). 
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